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Abstract 



Crack propagations in quasi-static fracture are studied theoretically. The 
Griffith theory is applied to discuss a crack extension condition and motion of 
crack tips in straight propagations. Stability of the straight propagations is 
investigated based on the simple assumption that a curvature near the crack 
tip is determined by a singular shear stress. It is shown that straight propa- 
gations become unstable under ceatain conditions. Combined with boundary 
effects, this instability leads to an oscillatory propagation of cracks. The 

critical conditions are calculated quantitatively. 
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Fracture has been studied extensively since Griffith wrote a breakthrough paper in 1920 
but still provides us with many unsolved problems 0. Yuse and Sano have recently 
carried out a nice experiment reminiscent of directional solidification [Q. In their experiment, 
a thin glass plate seeded with a small crack is moved from a hot region to a cold one at a 
velocity v. It was observed that various forms of crack patterns developed depending on the 
temperature difference 5T and the velocity v. When the temperature difference is increased 
and reaches a critical value, the crack begins to extend along a line. At a larger temperature 
difference, a transition from the straight crack to a wavy crack occurs. In this letter, we 
present a theory in which the two critical temperature differences are calculated. 

An experimental configuration is displayed in fig. ||. Here, a coordinate system (x, y) 
is fixed on the glass plate, which is assumed to be sufficiently thin, to have boundaries at 
y = ±W/2, to be submerged in a cold bath in the region — oo < x < s, and to contact a 
heating element at x — > oo. We note that s is regarded as a variable specifying the time- 
dependent external condition and that the propagation velocity of cracks is close to ds/dt. 
This fact implies that the crack tip moves adiabatically with respect to the change of the 
external condition. We thus parametrize the time-evolution of the position of the crack tip 
by s as (g(s),p(s)). The problem we address is then to derive an equation of motion for 

Suppose that a crack is positioned on the center line of a glass plate. According to the 
Griffith theory , when the strain-energy release rate G exceeds the energy needed to create 
a new surface F, the crack begins to extend at a speed of on the order of propagation of 
sound. For a fixed temperature profile, however, the crack would be arrested immediately 
at a position where G equals to F, because elastic forces arising from non- uniform thermal 
expansions are localized over a region where the temperature gradient is large. Then, if 
a temperature profile were translated infinitesimally along the plate, the crack tip would 
move because G > F, and then stop again at G = F after infinitesimal extension. Since the 
translation velocity v is much smaller than the sound velocity in glass, the crack propagates 
so as to maintain the equality G = T . Such a fracture is called 'equilibrium fracture' or 
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'quasi-static fracture'. This should be distinguished from an 'unstable fracture' or 'dynamic 
fracture' in which the energy release rate G becomes larger as the crack extends. In quasi- 
static fracture, the equation G = T gives not only a condition for the extension of a crack, 
but also the position of a tip of the straight crack. 

In order to embody this idea, we focus on infinitely narrow cracks. Then, stresses aij 
diverge at the crack tip. Using a polar coordinate system (r, 6'), we can expand cry in y/r 
around the crack tip. In particular, the singular part of the stress is characterized by stress 
intensity factors Ki and Ku related with basic 'modes' of crack-surface displacement called 
mode / (opening mode) and mode // (sliding mode), and its form is expressed by 



where t}i^{9) and Ylf^{9) are universal functions 0. Since the shear stress vanishes at 
the center of glass plates, the equality Kjj = holds for the straight propagation along 
the center line. Then, as shown by Irwin, the energy release rate G is related to the stress 
intensity factor Kj as G = K]/E, where E is a Young modulus |^. By introducing a critical 
stress intensity factor Kj such that Kj = VTE, the equation G = F is reduced to 

Kj = KI, (2) 

where Ki is determined by the crack pattern and temperature profile, and K'j is intrinsic 
to the material in question. Noting that Ki is a function of q{s) — s due to translational 
symmetry in the x-direction and that Kj does not depend on s because of (|^), we obtain 

q{s) = s + qo, (3) 

where go is a constant given by (^. A crack can not extend when a solution of (^ does 
not exist. The expression of Kj for a straight crack pattern will be given below, but before 
giving this derivation, we will discuss an instability of the straight propagation in the next 
two paragraphs. 

To this point, the system has been assumed to be under pure mode I loading. However, 
once the crack is inclined or away from the center, a shear mode is superposed onto mode 



/, and the crack is deflected from straiglit line geometry. We assume here that the crack 
extension condition remains unchanged when the crack extends smoothly. Then, in 
order to describe the dynamics of tips of curved cracks, a relation between shear stresses 
and curvatures at the crack tips should be considered. Our hypothesis is as follows: a 
curvature at a tip is determined by only and the crack tip tends toward the orientation 
of minimum shear loading. Since the curvature is well approximated by (Pp/ds^ for a slightly 
curved crack, dimensional analysis leads to a simple expression: 



where the non-dimensional function ^{z) is odd in z and negative for 2; > 0, and a new 
parameter rg with a length dimension has been introduced. A similar equation was consid- 
ered by R.V. Gol'dstein and R.L. Salganic 0. They identified tq as the length of the end 
region of cracks and derived the equation Kn = by taking the limit rQdp^/ds^ 0. Our 
investigation is not limited to this limit. In fact, we will now derive the dependence of the 
crack oscillation on Tq. 

Since Kjj{s) depends on the crack pattern which is identical to the history of the crack 
tip, {q{s') , p{s')) , (—00 < s' < s), (^) is regarded as a non-Markov dynamical system. We 
thus need to consider perturbations with infinite degrees of freedom so as to investigate 
stability of the straight propagation of cracks. However, if instabilities are caused by per- 
turbations with a longer curvature radius than the size of a most relevant region to Kn 
near the tip, we are allowed to employ a Markov approximation. Perturbations are then 
restricted to a linear combination of 'tilt' and 'translation' of the crack, and Kjj is written 
as 



where the first and second term of the right-hand side correspond to the stress intensity 




(4) 



(5) 



factors of shear stresses coming from 'tilt' and 'translation', respectively. Substituting 
into we obtain a linearized equation around the straight propagation: 
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where we have used a normahzation of $ such that $'(0) = —1. The oscillatory instabihty 
observed in experiments then can be explained by checking the change of the sign of K^^^ 
and the positivity of Kfj. The critical point is given by Kfj — and therefore does not 
depend on the unknown parameter tq, while the wavelength of oscillation A is proportional 
to 



^ (7) 
(0)- V) 



Our theoretical result about the critical point can be compared to that in experiments and 
the length tq can be estimated from experimental data. 

Now, we derive expressions for fC/, Kf^ and K^^. Since glass plates arc assumed to 
be sufficiently thin, we consider two-dimensional thermo-elastic systems under plane stress 
conditions. The equation of the force balance djCij — can be satisfied automatically 
when the stresses are written in terms of the Airy stress function (p such that 

Then, using Hook's law and the compatibility condition of a set of strains, we obtain a basic 
equation in our system: 

AV = -EaAT. (9) 

Here, o; is a thermal expansion rate and the temperature profile T{x, t) along the glass plate 
is derived under the assumptions that T{x, t) equals to Tc in the cold bath x < s, and that 
T{x, t) obeys a diffusion equation dtT = DAT under a boundary condition T{x, t) -^T^ + ST 
for X ^ oo, where L> is a thermal diffusion constant. A steady moving solution at a velocity 
V is obtained in the form 

T{x,t)^Tc + e{x- s)5T{l-ex^{-{x- s)/dQ)), (10) 
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where we have introduced a thermal diffusive length: do = D/v. Stress intensity factors 
are given from solutions of under boundary conditions J2j '^■ij^j = 0; where rij is a 
normal unit vector to boundaries. We solve the problem by dividing it into two steps. First, 
we calculate the stress dij in glass plates without cracks, where boundary conditions at 
y = ±W/2 and |x| oo are assumed as dyy = a^y = at y = ±W/2 and a^x = CTxy for 
oo. The following expressions (with cr° (x) denoting aij{x,y = 0)) are then obtained: 



\x 



aL(x) = -Ea / dke''^f{k)Uk), (11) 

J —oo 

and 

dke''''T{k){l-U{k)), (12) 

-oo 

where T{k) is a Fourier mode of the temperature field given by (p!0| ) and 

_ cosh{kW/2) ± 2 smh{kW/2) 
^^^^^ ~ kW + sinh(A;iy) " ^^^^ 

Then, noting that existence of a crack alters only boundary conditions, we find that the 
stresses on the glass plate with a crack are given by aij + 0"*^, where a,*- is a solution of the 
homogeneous equation of (P) under fictions external forces introduced so as to satisfy the 
boundary conditions along the crack. Since a singularity appears in the stresses 0"^*-, stress 
intensity factors Ki and Kn for a straight crack are expressed by the fictions forces such 
that [| 

= \ - I dr^, and Kn = J- f dr^, (14) 
\ n Jo Jr \ 71 Jo \/r 



where dee and are correspond to the fictions forces contributing to Kj and Kjj respectively, 
and r is measured along the crack from the tip. This formula, strictly speaking, should be 
applicable to cracks extending semi- infinitely in an infinitely large space. We expect however 
that this approximation is still effective for sufficiently large W/do. By carrying out a partial 
integral, we obtain the expression of Kj for a straight crack on the line y = 0: 

Ki = t\^ r dx{dlal){q-xf/\ (15) 
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Also, noting a transformation (Jre = o'xy — {(Jxx — cryy)dp/ds, we can obtain Ku for cracks 
inclined slightly and displaced from the center: 

K,. = /' + ^ r (16) 

V TT j-oo ^/q ~ X as J-oo y/q — x 

where the stresses have x-dependence in the form aij{x,y{x)), y{x) = p + dp/ds ■ {p — x). 
We further expand the right-hand side of this equation in p and dp/ds. The expansion 
coefficients K^^ and Kf'^ are then derived as 

4? = -2^/^r rfx(9,VL)v/^, (17) 

V TT J-oo 

and 

Kf^ = Kj- 2./? /' rfx(9^aL) ■ (q - xfl\ (18) 

V TT J-oo 

where we have used the equation of force balance: dya^y = —dx^xx- From (|lT]),(|T2D, (^5|), 
(0) and (|18|), we obtain fsT/, i^i?, and x'ff . 

We now introduce two control parameters defined by /i = W/d^ and i? = E'a^Tl^^/^/i^j. 
One may easily see that /i and R correspond to non-dimensional velocity and temperature 
difference, respectively. We further define ^ as a ratio of a crack length out of the cold 
bath to a thermal diffusive length: ^ = qo/dQ. Then, Kj, Kj^j and Kj^j are expressed 
as Ki = K'^RFiifi,^), Kf^ = K^RfP {fi,^)/W , and K^^j' = K'^RFl]\fi, where Fj, 
F^]^ are non-dimensional functions calculated numerically. When we compute these 
numerical values, we fix the parameter value /i = 10 as representaion of ideal experimental 
conditions under which the thermal diffusive length is much larger than a thickness of glass 
plates. 

For sufficiently small R, there is no solution ^ satisfying RFj{fi, = 1 which is equivalent 
to Kj = Kj. In this case, the crack can not extend. For R > R^^\ RFj(fi,^) = 1 has a 
solution ^ as is shown in fig. ^ Then, when the temperature difference is increased, the 
crack length increases, and Fj]^ changes its sign at R^^^. The straight propagation is stable 
for R'^^ < R < R'^^\ but the shear stress arising from tilt makes the crack more inclined 



for R > R^^^ ■ We also confirmed the positivity of F^^^ . Therefore, we can say that a 
Hopf bifurcation occurs at Rf"* with a critical wavelength Ac- Here, R^p , Rf^ and Ac were 
computed as R^p = 5.6, R'^'^ = 9.9, and Ac = 4:.5y/roW. Similarly, for an arbitrary fx, 
assuming that ^/(/i, ^) takes a maximum value at ^m{l^) and that Fj]\ij,,^) becomes zero 
at ^o{^), we obtain 

= ^ Fj{fx,Uf-^))-\ (19) 

and 

^ \ f(o)^ c ( ^n V^o^- (20) 
^ ^11 if^,W)) 

Let us discuss asymptotic forms of these equations for fi oo. First, from (0) and (|12|), 
we see that the stress cr° is proportional to fi in this limit. Then, ([l5|) , ([T7|) , (pTsD , (p!9D , and 
(pop lead to asymptotic forms for /i — > oo: 

^ i?(2) ^ ^-1, and Ac ~ (21) 

By using dimensional quantities, the asymptotic forms are rewritten as 5Tc ~ and 
Ac ~ f° for V oo, and 5Tc ~ and Ac ~ VW for W oo. The latter relation 

implies that the oscillatory instability is caused by effects of the boundaries. In fact, if we 
assume the periodic boundary conditions at y = ±W/2, Kjf = is obtained due to the 
translational symmetry in the y-direction and the crack then propagates at a finite angle 
to the line y = when K^j < 0. Similar behavior was observed in experiments using glass 
cylinders instead of glass plates. 

In the experiment done by Yuse and Sano 0] , the thickness of the glass plates is greater 
than the thermal diffusive length in the range f > 5 [mm/s]. Therefore, their experimental 
system may be regarded as a two-dimensional one only in the range f <^ 5 [mm/s]. The 
relation (|21| ) seems to be consistent with experimental results in the range 0.5 [mm/s] < 
f ^ 5 [mm/s]. In order to see correspondence with experiments more quantitatively, we use 
the following values of the material constants typical for glass plates: a = 8 x 10~^ [K~^], 
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E = 7 X 10^° [Pa], r = 8 [J ■ m-2], D = 5 x lO'^ [cmVs]. Then, for glass plates with 
W = 2.4 [cm], n = 10 corresponds to f = 0.2 [mm/s], and the critical temperatures 5TJ^) 
and 6TP are computed as 6T^^^ = 50 [K] and 6TP = 80 [K]. In the experiment @], 
was found to be somewhat less than our result, while Tp'' was approximately twice of ours. 
We believe that this discrepancy may be due to error in estimating F. It is not yet known 
how to accurately measure this value. Further investigations will be necessary. 

Finally, we briefly discuss the parameter ro introduced above. From the experimental 
data that Ac = 3 [mm] for W = 1 [cm], we obtain tq = 0.05 [mm]. As we expected above, 
this value is much smaller than curvatures of oscillating cracks. At this time, however, we 
are not able to satisfactly estimate the value of tq using physical arguments. By clarifying 
what occurs on this scale, we hope to reach a deeper understanding of fracture phenomena. 

In summary, we have proposed a theory for the transitions between the crack patterns in 
quasi-static fracture. Based on our theory, the critical points for yU = 10 have been calculated 
quantitatively and the scaling relations for yU — > oo have been derived. Recently, M. Marder 
attempted to explain the oscillatory instability [^ relying on the theory of Cotterell and 
Rice []TU|. In addition, numerical simulations of discrete models were carried out []TT|. We 
have not yet been able to work out the correspondence between these works and ours. In 
order to obtain a true understanding of this phenomena, we need to synthesize theories, 
numerical simulations and experiments. 
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FIGURES 

FIG. 1. Schematic figure of an experimental configuration. The glass plate is submerged in the 

cold bath represented by the dotted region. 



FIG. 2. ^ versus R (sohd line) and K^^/Kj versus R (dotted line) . Vertical axes of solid and 
dotted lines are graduated on the left and right, respectively. These graphs starts from R^\ and 
Kf^/Kj is negative for R > R^ . Here, R^^^ = 5.6, and R^ = 9.9. 
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